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In this letter, we propose a model of inflationary cosmology with a bounce preceded and study 
its primordial curvature perturbations. Our model gives rise to a primordial power spectrum with a 
feature of oscillation on large scales compared with the nearly scale-invariant spectrum generated by 
the traditional slow rolling inflation model. We will show this effect changes the Cosmic Microwave 
Background (CMB) temperature power spectrum and the Large Scale Structure (LSS) matter power 
spectrum. And further with a detailed simulation we will point out this signal is detectable to the 
forthcoming observations, such as PLANCK and LAMOST. 



Inflation, as a description of the very early universe, 
has successfully resolved some problems existing in hot 
Big Bang cosmology, such as flatness, horizon, monopole 
problem and so on However, this scenario is puz- 
zled by the initial singularity [2|. One possible approach 
to this disaster is to introduce a bounce before the in- 
flationary expansion, which requires the hot Big Bang 
expansion be preceded by a contracting period js], 0, [B[ . 
If this happens, one significant question would be pro- 
posed: what does bouncing cosmology tell us for obser- 
vations? Or, is it detectable for some primordial relics 
from contracting phase to be imprinted on observations? 
To answer this question, we need to study the evolution 
of primordial gravitational perturbations seeded before 
the bounce. 

In this letter, we propose a nonsingular inflationary 
model, and study signatures of its primordial fluctuations 
on CMB temperature power spectrum and LSS matter 
power spectrum. We find an interesting oscillation signa- 
ture existing on large scales of the scale-invariant spec- 
trum and by a detailed simulation we will show this new 
effect could be detected by the forthcoming astronomical 
observations, such as PLANCK and LAMOST. 

As in inflation theory, our model can be described in 
terms of scalar fields which minimally couple to the four 
dimensional Einstein's gravity. Explicitly it consists of 
two scalar fields <f> and if) with the lagrangian given by: 
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in a spatially flat Friedmann-Robertson- Walker (FRW) 
universe. Here the essential component is the scalar field 
ijj. It plays a crucial role in giving a bouncing solution 
smoothly. Without it, the model in Eq. (1) will be sim- 
ilar to the traditional inflation with a single scalar field, 
which as we know suffers from the problem of the ini- 
tial singularity. In this model the potential is only the 
function of the field <f> and of Coleman- Weinberg form Q : 



(2) 



which takes its maximum value Aw 4 /16 at 4> = and van- 
ishes at the minima when <b — ±v. Therefore, the scalar 



field ip merely affects the evolution around the bounce 
but decays out quickly when away from it. 

In order to discuss the perturbations explicitly, we first 
see how the background universe evolves. In this model 
a contracting universe can be driven to reach a mini- 
mal size during which the universe evolves like a matter- 
dominant one, and then a quasi-exponential expansion 
is following, and so is able to explain the problems ap- 
peared in standard Big Bang cosmology. The process to 
link the contraction and expansion is a smooth bounce, 
and the evolution of the hubble parameter can be treated 
as a linear function of the cosmic time approximately. 

We take the initial condition for the background as 
that (j) stays at one vacuum like — v when the universe is 
contracting and ij> is small enough which can be ignored 
on background evolution. In this phase, the field cf> os- 
cillates around — v making the equation-of-state (EoS) of 
the universe oscillate about w = 0, and so the average 
state being similar to a matter-dominated one. Thus we 
have the useful expressions of background evolution 



a ~(- v ) 2 , n = - , \4>\^ v - 
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where TL = a' /a is the comoving hubble parameter and 
the prime denotes the derivative with respect to the co- 
moving time. Another useful relation is given by 



2TLTL' - TL" 
2(H 2 - TL') 



(4) 



which will be used to calculate the metric perturbations. 

Since the universe is contracting, the amplitude of the 
field 0's oscillation gets larger and larger, while the con- 
tribution of the field ip grows rapidly. When the field 
reaches the plateau, the bounce happens at the mo- 
ment During the bounce, we take the parametriza- 
tion H(t) = a(t — tg) around the bounce point tg, 
and the coefficient a is a positive constant determined 
from numerical calculations. In the bouncing phase, 
the kinetic term of ip reaches the maximal value and 
from the equation of motion we deduce an expression 
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-3H when a is not very 
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large. Finally we have the approximate relations 
H a |(77 - VB ), r ~ 1^1 ~ e-Jf(™) 2 , (5) 

where we have defined y = 8aa^/ir. 

After the bounce, as the field 4> moves forward slowly 
along the plateau, the universe enters into an expanding 
phase at the moment ts+ and the EoS of the universe is 
approximately — 1. The universe expands with its scale 
factor growing almost exponentially. In this phase, we 
have the well-known relations for background evolution 

a ~ , H ~ Constant . (6) 

V 

Finally, when the field drops into the vacuum +v, it will 
oscillate again and the EoS of the universe will oscil- 
late around zero as it does before the bounce. To make 
the scenario explicit, we give a sketch description of our 
model in Fig[T] and present the numerical calculation of 
the background parameters in Fig [2] 




FIG. 1: Evolution of the universe. A sketch plot of the 
evolution of the universe in the model of Eq.fTJ. Before infla- 
tion, there is a bounce instead of the initial singularity. 

Now we study the linear perturbations of the model. 
Taking the longitudinal (conformal Newtonian) gauge, 
the metric perturbation is presented as follows: 

ds 2 = a 2 {r 1 )[(l + 2^)dr 1 2 -(l-2^)dx i dx l ] , (7) 

and the equation of motion of the gravitational potential 
is: 

$" + 2(H - + 2{H' - ^4-)$ - V 2 $ 

<p' qy 

= 8irG(2H + ^jWH , (8) 
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FIG. 2: Evolutions of <f>, w and H. A plot of the evolutions 
of the field <f>, the EoS w and the hubble parameter H in the 
model of Eq.(JT| where we take the Coleman- Weinberg poten- 
tial. In the numerical calculation we choose the parameter 
A = 8.0 x 1CT 14 , v = 0.82M pi , and the initial condition as: 
4> = -0.82M pi , <j> = 3.0 x 10 _:L0 M^, ip = -0.72M pi , tji = 
5.0 x W' 13 M^ where M pt = 1/VC . 



which can be derived from the basic perturbation equa- 
tions directly (we refer the complete derivation to Ref . 0] , 
and see e.g. Ref.[8j for a comprehensive survey of the cos- 
mological perturbation theory). As is pointed out previ- 
ously, the energy density of the field "0 is usually negligi- 
ble far away from the bounce, and hence we have ip' ~ 0. 
Near the bounce, ip becomes very important, but accord- 
ing to the analysis in Ref.[7j we have the approximation 
2H + <fi"/(j>' ~ and so the perturbation of ip decouples 
from Eq. |8]). Therefore, we will neglect the r.h.s. of 
Eq. ©, and just focus on the adiabatic fluctuations in 
the following which can be determined by a single scalar 
field (j). 

Now we follow one Fourier mode of the perturbation, 
labelled by its comoving wave number k, and find that 
there are two paths for perturbations. The evolution of 
perturbations is sketched in FigJ3] Initially all the pertur- 
bations stay inside the horizon in the far past. Since the 
hubble radius shrinks, those modes with small comoving 
wave number exit the horizon while the large k scales still 
keep inside. When the bounce takes place, all the pertur- 
bations will enter the horizon because at that moment the 
hubble radius diverges. Since in our model the bounce is 
followed by an expanding phase, those Fourier modes will 
escape out if the efolds for the post-bounce slow-rolling 
period is large enough. After that, these modes will re- 
enter the horizon at late times after the slow-rolling phase 
has finished. In the following calculations, we will focus 
on large k region and see whether the large k modes are 
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FIG. 3: Evolution of perturbations. A sketch plot of 
the evolution of perturbations with different comoving wave 
number k in our model. 



able to perform scale-invariant spectra and give more in- 
formation on the CMB observations. 

For the contracting phase before the bounce, the equa- 
tion of motion in momentum space can be solved ex- 
plicitly. We take the Bunch-Davies vacuum as the ini- 
tial condition <j>^ ~ \<p\e^ tkv when the perturbation 



are deeply inside the horizon. 
3 we obtain 



Since during this period 
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where the subscript i represents the initial time. Substi- 
tuting the Eq. © into Eq. ((8|) and solving it, we have 
the solution to the perturbation in the bouncing phase 



\v(v-vb) 



x I C fe cos[fc(?7 - r) B )] + D k sin[fc(?7 - r} B )] M10) 



Moreover, for the nearly de-Sitter expanding phase, we 
obtain 



k v F k J_ v {k{r,- 



v E k J v {k^- 
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where 7 ~ 1/2, v ~ 1/2 and fjB+ = Vb+ + 1/"Hb+- 

Having obtained the solutions of the perturbation 
equations in different phases, it is necessary to know 
the matching relations among these solutions and deter- 
mine the coefficients C k , D k , E k and F k . This depends 
on whether the curvature perturbation on a uniform co- 
moving hypersurface or the gravitational potential passes 



through the bounce regularly [§] (see also [l(| for a re- 
cent study). For a nonsingular bounce scenario such as 
what we considered, the continuity of background evolu- 
tion implies that both $ and $' are able to pass through 
the bounce smoothly. By matching $ in Eqs. © and 
(JTUJ) on the surface tjb-, and that in Eqs. ([TO]) and (fTTj) 
in sub-hubble region on the surface t\b+ as well as their 
comoving time derivatives, all of the coefficients can be 
determined. However, since E k represents a decaying 
mode when escape outside the horizon, we neglect it and 
finally obtain the dominant mode 



F k 




sm[k(r, B - - fj B +)] H12) 
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By comparing the coefficient of Eq. (|T2|) and the initial 
form of <I> in Eq. |9]), one obviously notice that the sub- 
hubble form of the metric perturbation has obtained an 
oscillation term when the universe undergoes a bounce. 
Note that, another important quantity is the curvature 

7 (*' + 
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perturbation in comoving coordinate £ = 3>+ 
7i$), and when the expansion is nearly de-sitter like, 
there is a simple relation between these two quantity £ ~ 
$/e with the slow roll parameter e = —H/H 2 . Therefore, 
we eventually have the primordial power spectrum for the 
curvature perturbation 
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where we have assumed t}b+ ~Vb- to be small compared 
with I ^ |. Obviously, the first term provides a nearly 
scale-invariant spectrum which is consistent with current 
cosmological observations. However, the second term ap- 
parently shows that there is a wiggle on the spectrum, 
due to the modification brought by a bounce. 

For a numerical estimate of TLb+ and TIb-, we nor- 
malize the current scale factor ao = 1 and choose the 
current hubble parameter Hq = 72 km s _1 Mpc -1 = 
1.536 x 10~ 42 GeV, the hubble parameter during infla- 
tion Hi ~ 1.68 x 10 12 GeV, and the e-folds for inflation 
N ~ 60. Therefore, our model predicts that TCb+ — 
2 x 10" 4 Mpc _1 and Ub- ^ -1.6 x lO^Mpc" 1 . Based 
on the primordial spectrum in Eq. (|13p. in Fig[4]we illus- 
trate the CMB temperature power spectrum and present 
LSS matter power spectrum. For comparison we have 
also considered the standard case where TLb- is taken to 
be zero. In the numerical calculations we use the publicly 
available Markov Chain Monte Carlo (MCMC) package 
CosmoMC [ri| and take the basic cosmological parame- 
ters as given below: 



{n b h 2 ,n c h 2 ,T,H ,A s ) 

(0.022, 0.115, 0.088, 72, 2.3 x 10~ 9 ) 



(14) 
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and the pivot scale fc* = 0.05 Mpc -1 . One can see from 
Fig|4] that our model leads to an obvious fc-dependent 
oscillation signature in the power spectrum, especially at 
large scales. 
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it is argued that fluctuations before and after the bounce 
are largely independent (yet see Ref.[l9| for some 
criticisms). In this letter we propose a concrete cosmo- 
logical model with inflation preceded by a bounce, and 
by investigating in detail the perturbations we show some 
imprints of the bounce are detectable to the forthcoming 
CMB and LSS observations. 



Acknowledgments 

We have performed our numerical analysis in the 
Shanghai Supercomputer Center (SSC). We thank 
Robert Brandenberger, Mingzhe Li and Paul Steinhardt 
for useful comments on the manuscript. This work is 
supported in part by National Science Foundation of 
China under Grant No. 10533010, and the 973 program 
No.2007CB815401, and by the Key Grant Project of Chi- 
nese Ministry of Education (No. 305001). 



FIG. 4: Observational effects by a bounce. The effect 
on the CMB temperature power spectrum and LSS matter 
power spectrum by our model H.B+ — 2 x 10 _4 Mpc _1 and 
TLb- ~ —1.6 x 10~ 4 Mpc _1 . The dots and error bars are 
WMAP5 and SDSS data. 

To test our model we firstly consider the current astro- 
nomical observations from WMAP5 12] and SDSS [H, 
due to the large uncertainties at large scales as shown 
in FigU] by the blue dots and the error bars, we find 
that the oscillating spectrum of our model is consistent 
with the data. Thus we consider the forthcoming mea- 
surements PLANCK [3] and LAMOST with higher 
precision. We simulate the CMB TT, TE and EE power 
spectra with the sensitivity of PLANCK and the LSS lin- 
ear matter power spectrum with the sensitivity of LAM- 
OST and find these measurements will be sensitive to 



TLb- 



-7.0 x 10 5 Mpc which is smaller than our 



predicted value and makes our model testable. If this 
signal would be detected, it will act as a smoking gun to 
the bouncing cosmology. 

Physics of bouncing cosmology, since it happens in ex- 
tremely high energy regime, is hardly to be found by 
experiments directly. So it is a debate whether a bounce 
has taken place or not. To find the evidences of a bounce, 
we need to know what can a bounce leave for observa- 
tions. This question is still discussed drastically in the 
literature, and one potential clue is to study the primor- 
dial curvature fluctuations. In the context of the Pre- 
Big-Bang scenario |3j and in the cyclic/Ekpyrotic cos- 
mology [4j , the resulting curvature perturbation strongly 
depends on the physics at the epoch of thermalization, 
and thus an uncertainty of a thermalized surface is in- 
volved lj| 17 1. In the frame of loop quantum cosmology, 
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